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HIGHLIGHTS 


•  Improved  agglomerate  model  is  used  to  study  the  PEM  fuel  cell  catalyst  layer. 

•  Catalyst  and  ionomer  loadings  varied  in  through-thickness  and  in-plane  directions. 

•  Higher  catalyst,  ionomer  loading  at  membrane/CL  interface  improves  performance. 

•  Bidirectionally-graded  CLs  are  investigated  for  the  first  time. 

•  Bidirectional  grading  improves  performance  in  ohmic  and  mass  transport  regimes. 
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The  catalyst  layer  (CL)  of  the  polymer  electrolyte  membrane  (PEM)  fuel  cell  must  be  modeled  accurately 
in  order  to  resolve  the  effects  of  complex  interactions  between  charge  and  mass  transport  on  the  fuel 
cell's  electrochemical  reactions.  In  previous  work,  we  developed  an  agglomerate  model  [1]  which 
correctly  accounts  for  variations  in  the  agglomerate  surface  area  as  the  CL  constituents  are  varied  to 
provide  a  better  estimate  of  diffusion  losses.  Here,  this  improved  agglomerate  model  is  employed  to 
investigate  a  PEM  fuel  cell  catalyst  layer  with  a  functionally-graded  composition.  We  present  results  for 
varying  catalyst  and  ionomer  loadings  in  both  the  through-thickness  and  in-plane  directions.  In  agree¬ 
ment  with  experimental  observations,  we  find  that  a  higher  catalyst  and/or  ionomer  loading  at  the 
membrane/CL  interface  improves  performance  especially  in  the  ohmic  loss  regime.  Similarly,  improved 
performance  is  observed  for  higher  catalyst  and/or  ionomer  loadings  under  the  channel  in  the  mass 
transport  loss  regime.  In  addition,  we  investigated  bidirectionally  graded  CLs  for  the  first  time.  It  is 
observed  that  higher  performance  can  be  obtained  with  bidirectionally  graded  CLs  in  both  ohmic  and 
mass  transport  loss  regimes. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Over  the  past  two  decades,  many  important  developments  have 
brought  PEM  fuel  cells  closer  to  commercialization.  Although 
substantial  experimental  and  computational  research  has  been 
conducted  on  all  aspects  of  PEM  fuel  cells,  the  sluggish  reactions  in 
the  cathode  catalyst  layer  (CL)  remain  a  major  limiter  of  perfor¬ 
mance.  Because  all  of  the  fuel  cell's  electrochemical  reactions  and 
significant  charge  and  mass  transport  occur  within  the  catalyst 
layer,  the  overall  accuracy  of  fuel  cell  computational  models  de¬ 
pends  critically  on  an  accurate  representation  of  the  catalyst  layer. 
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Therefore,  modeling  and  simulation  of  the  CL  continues  to  attract 
significant  research  effort. 

The  CL  is  a  10-30  pm  thick  composite  structure  containing 
platinum  (Pt),  carbon  support,  ionomer  electrolyte,  gas  pores,  and 
liquid  water.  Within  this  structure,  reactions  are  catalyzed  by  Pt 
while  the  carbon  support  and  ionomer  phase  conduct  electrons  and 
protons,  respectively,  and  the  pores  create  pathways  for  reactant 
gas  access  and  product  removal.  Consequently,  the  CL  activity 
consists  of  multiple  competing  transport  processes  which  require 
an  optimum  loading  of  catalyst,  carbon  support,  and  ionomer 
distributed  within  a  layer  of  adequate  porosity.  For  example, 
whereas  a  higher  Pt  loading  would  improve  the  reaction  rate,  it 
would  adversely  impact  reactant  transport  due  to  decreased 
porosity.  Similarly,  a  higher  volume  fraction  of  ionomer  would 
improve  proton  transport,  but  also  at  the  cost  of  reduced  reactant 
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transport  due  to  a  lower  porosity.  In  this  context,  several  experi¬ 
mental  studies  [2-5]  have  been  performed  to  optimize  the  CL 
composition  to  improve  the  overall  PEM  fuel  cell  efficiency. 

Uchida  et  al.  [2]  experimentally  studied  the  effects  of  ionomer 
loading  and  reported  an  optimum  value  for  ionomer  loading  that 
gave  better  performance  over  the  entire  operating  regime.  Simi¬ 
larly,  Antolini  et  al.  [3]  reported  an  optimum  Nafion  loading  which 
minimizes  activation  and  ohmic  losses.  Paganin  et  al.  [4]  examined 
both  platinum-to-carbon  weight  ratio  (%Pt|C)  and  Nafion  loading 
and  also  reported  optimum  values  for  these  quantities.  Passalacqua 
et  al.  [5]  unified  the  optimum  Nafion  loading  results  of  [2-4]  by 
defining  a  parameter  called  NFP  (Nafion  percentage).  They  also 
performed  experiments  and  reported  an  optimum  NFP  value  of  33% 
which  agrees  with  the  optimum  NFP  loading  of  about  33-40% 
given  by  Refs.  [2-4].  Qi  and  Kaufman  6]  studied  the  effects  of  Pt 
loading  and  %Pt|C  for  a  CL  with  33%  NFP  which  was  the  optimal 
Nafion  loading  according  to  their  experimental  data.  They  reported 
optimum  Pt  loadings  for  20%  and  40%  Pt|C  ratios. 

The  literature  contains  many  other  experimental  studies  that 
focus  on  the  optimization  of  the  CL  composition.  However,  most  of 
these  consider  a  spatially-uniform  CL  composition.  There  are  rela¬ 
tively  few  experimental  attempts  [7-11]  to  investigate 
functionally-graded  CL  compositions.  Ticianelli  et  al,  Antoine  et  al, 
[7-9  studied  the  effect  of  a  graded  catalyst  loading  through  the 
thickness  of  the  CL.  Ticinalli  et  al.  [7,8]  reported  performance 
improvement  when  the  catalyst  concentration  was  increased  at  the 
membrane/CL  interface.  Antoine  et  al.  9]  investigated  porous  and 
non-porous  active  layers.  For  the  porous  case,  they  observed  per¬ 
formance  improvement  when  the  catalyst  loading  was  increased 
adjacent  to  the  membrane.  On  the  other  hand,  they  reported  better 
performance  when  the  catalyst  loading  was  higher  at  the  gas 
diffusion  layer  (GDL)/CL  interface  for  the  non-porous  electrode. 
Wilkinson  and  Pierre  10]  studied  graded  catalyst  loading  along  the 
length  of  the  gas  channel  and  reported  enhanced  performance 
when  the  catalyst  loading  was  higher  towards  the  inlet.  In  addition, 
Xie  et  al.  [11]  examined  the  effect  of  an  NFP  gradient  through  the 
thickness  of  the  CL  and  reported  improved  performance  for  a 
higher  NFP  towards  the  membrane/CL  interface. 

Determining  an  optimal  CL  composition  requires  a  good  un¬ 
derstanding  of  the  physical  phenomena  occurring  within  the  CL. 
Numerical  models  serve  as  a  valuable  tool  to  provide  detailed 
insight  that  might  be  difficult  to  observe  experimentally.  The 
literature  contains  several  numerical  studies  [12-17]  that  have 
investigated  functionally-graded  CL  compositions.  Wang  et  al.  [12] 
developed  a  ID  cathode  model  based  on  percolation  theory  that 
characterizes  effective  CL  properties  and  the  active  reaction  area. 
Employing  a  macro-homogeneous  model  for  the  CL,  they  investi¬ 
gated  the  effects  of  graded  NFP  through  the  CL  thickness.  In 
agreement  with  experimental  results  [11],  they  reported  better 
performance  for  high  NFP  close  to  the  membrane/CL  interface. 
However,  their  use  of  the  macro-homogeneous  approach  to  model 
the  CL  prevented  them  from  predicting  the  effects  of  the  graded  CL 
in  the  mass  transport  regime.  Using  the  classical  agglomerate 
approach,  Sun  et  al.  13]  reported  improved  performance  for  higher 
Nafion  loading  towards  the  membrane  with  their  2D  cathode 
model.  Song  et  al.  14]  optimized  the  Pt  and  NFP  distributions 
through  the  thickness  of  the  CL  with  their  ID  cathode  model  which 
employs  a  macro-homogeneous  approach  by  considering  the  effect 
of  liquid  water  saturation.  They  obtained  the  optimum  NFP  distri¬ 
bution  as  a  linearly  increasing  function  from  the  membrane  to  the 
GDL.  For  the  optimum  Pt  loading  distribution,  they  obtained  a 
convex  increasing  function  from  the  membrane  to  the  GDL. 

Schwarz  and  Djilali  15]  developed  a  multiple  thin-film 
agglomerate  model  to  account  for  the  effect  of  liquid  water  in  the 
CL.  With  their  improved  agglomerate  approach  they  investigated 


graded  catalyst  loading  cases  in  three  separate  directions  with  a  3D 
fuel  cell  model.  They  reported  improved  performance  when  the 
catalyst  loading  was  higher  towards  the  channel  inlet,  and  when  it 
was  higher  towards  the  channel  in  the  lateral  in-plane  direction. 
Srinivasarao  et  al.  [16]  investigated  the  graded  CL  composition  in 
the  through-thickness  direction  by  considering  multiple  stacked 
CLs  within  their  2D  cathode  model.  They  used  the  multiple  thin- 
film  agglomerate  model  as  in  Ref.  15]  and  investigated  graded 
catalyst,  NFP  and  Pt|C  distributions  separately.  In  addition, 
Mukherjee  and  Wang  17]  studied  a  bi-layer  CL  with  direct  nu¬ 
merical  simulation.  In  agreement  with  the  previous  studies  they 
also  reported  increased  performance  when  the  Nafion  loading  was 
higher  towards  the  membrane/CL  interface. 

It  should  be  noted  that  the  fidelity  of  PEM  fuel  cell  simulations 
depends  critically  on  the  employed  CL  modeling  approach.  The 
agglomerate  approach  is  the  most  comprehensive  among  the 
various  CL  models  which  are  reviewed  in  detail  in  Ref.  18]. 
Recently,  the  agglomerate  model  was  modified  by  employing 
discrete  catalyst  particles  [18,19]  in  order  to  more  accurately  predict 
the  effect  of  catalyst  loading.  In  addition,  in  Ref.  1  ,  we  presented 
an  improved  agglomerate  model  incorporating  a  sphere-packing 
approach  to  physically  model  the  effect  of  porosity  on  the  effec¬ 
tive  surface  area  of  agglomerates  (which  was  ignored  in  the  pre¬ 
vious  approaches)  as  this  is  a  crucial  parameter  to  determine 
diffusion  losses.  It  was  shown  in  Ref.  [1  that  the  new  agglomerate 
model's  predictions  were  physically  realistic  and  gave  good 
agreement  with  experimentally  observed  trends. 

In  this  study,  we  employed  the  model  in  Ref.  [1]  to  investigate 
the  effects  of  graded  CL  compositions.  In-plane,  through-thickness, 
and  bidirectionally-graded  CLs  were  studied  as  depicted  schemat¬ 
ically  in  Fig.  1.  It  should  be  noted  that  the  effect  of  CL  porosity  on 
surface  area,  which  is  captured  in  Ref.  [1  ,  was  not  accounted  for  in 
previous  agglomerate  models  13,15,16  while  investigating  the 
graded  CL  composition.  Moreover,  none  of  the  previous  studies 
provided  guidelines  for  the  range  of  variation  in  loading  within  the 
graded  CL.  Most  importantly,  the  multidirectionally-graded  CL 
composition  as  shown  in  Fig.  lc)  has  not  been  studied  to  date. 

2.  Model  description 

The  complete  2D,  steady-state  cathode  model  is  explained  in 
detail  in  Ref.  1  .  The  model  accounts  for  the  transport  of  chemical 
species,  electrons,  and  ions  in  the  GDL  and  the  CL  under  the 
assumption  of  an  isothermal,  isobaric,  and  steady-state  process. 
The  computational  domain  encompasses  half  of  the  channel  and 
half  of  the  land  area  as  shown  in  Fig.  2.  The  associated  geometric 
parameters  are  listed  in  Table  1. 

The  governing  equations  for  species  and  charge  transport  are 

[1]: 


-R> 

k  k 

(1) 

(<4ffV0m)  =  V-i 

(2) 

(<rfV0s)  =  -Vi 

(3) 

where  p  is  the  mixture  density,  (Oj  is  the  mass  fraction,  DJkff  is  the 
effective  diffusivity,  M,  is  the  molecular  weight,  Rj  is  the  source  term 
calculated  from  the  electrochemical  reactions,  is  the  effective 
ionic  conductivity,  0m  is  the  ionic  potential,  aff  is  the  electronic 
conductivity,  0S  is  the  electric  potential,  and  V  i  is  the  current 
generation  due  to  the  electrochemical  activity  in  the  CL  1  . 
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Fig.  1.  Schematic  representation  of  graded  CL  composition  in  various  directions. 


Fig.  2.  2D  computational  domain  with  half-cell  section  view  [1]. 


The  CL  is  modeled  by  employing  the  spherical  agglomerate 
approach.  As  described  in  Ref.  1  ]  the  sphere-packing  approach  was 
used  to  provide  as  accurate  estimate  of  the  effective  surface  area  as 
the  agglomerates  begin  to  overlap  with  decreasing  porosity.  The 
classical  agglomerate  equations  were  rearranged  with  the 
assumption  that  the  reactant  concentration  at  the  agglomerate 
surface  does  not  change  between  a  single  and  an  overlapping 
agglomerate.  Then  the  reaction  rate  is  calculated  based  on  the 
reactant  flux  at  the  agglomerate  surface  using  the  following 
equations  1]: 


surface,  Cs  is  the  oxygen  concentration  at  the  agglomerate  core 
surface,  ragg  is  the  agglomerate  radius,  5  is  the  ionomer  film  thick¬ 
ness,  aagg  is  the  effective  surface  area,  Dagg  is  the  effective  diffusivity 
in  the  agglomerate  core,  Th  is  the  Thiele  modulus,  and  H  is  Henry's 
constant. 

As  the  porosity  is  decreased,  adjacent  agglomerates  begin  to 
touch  and  then  overlap  which  reduces  the  available  agglomerate 
surface  area.  The  effective  surface  area  reduction  due  to  decreasing 
porosity  is  modeled  using  the  sphere-packing  approach  and 
calculated  as  [1]: 


Vi  =  I 


CL 


-4  FD 


(Cp  -  Cs) 
5 


aagg 


(4) 


a 


agg 


47rn(ragg  -T  (5)  for  ev  >  £Vjcr 


Cs 


_ Co _ 

1  +%s  (i^+s)  (coth(Th)Th  -  1^ 


where  F  is  the  Faraday  constant,  D  is  the  diffusivity  of  oxygen  in  the 
ionomer,  Co  is  the  oxygen  concentration  at  the  ionomer  film 


Table  1 

Geometric  parameters  for  the  2D  computational  domain. 


Channel  width 

wch 

1.5  mm 

Land  width 

W; 

1.5  mm 

GDL  thickness 

£gdl 

250  pm 

CL  thickness 

Cl 

15  pm 

aagg  — 


7T 


&mDp(]  ~  f) 


1 


mf 


for  ev  0  <  ev  <  £VjCr 


for  Sy  <  £p,0 


where  n  is  the  number  of  agglomerates  per  unit  volume  of  the  CL,  ev 
is  CL  porosity,  ev<cr  is  the  critical  porosity,  ev>0  is  the  limiting  porosity, 
Dp  is  the  agglomerate  diameter, /is  the  overlapping  ratio,  m  is  the 
number  of  contact  points,  and  am  is  the  packing  factor  1  .  Here  it  is 
important  to  note  that  the  catalyst  ink  composition  is  specified  in 
terms  of  the  weight  fractions  of  its  various  components  which  are 
expressed  in  terms  of  the  Pt  loading  to  calculate  the  porosity  as  in 
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Ref.  1].  The  am  and  m  values  employed  in  the  current  simulations 
correspond  to  the  rhombohedral  packing  scheme.  One  can  refer  to 
[1  for  additional  information. 

The  boundary  conditions  for  the  governing  equations  can  be 
written  by  referring  to  Fig.  2  as  follows.  For  species  transport 
(Maxwell-Stefan  equation),  the  inlet  mass  fractions  are  specified  at 
the  channel  and  the  no-flux  condition  is  assigned  to  all  the 
remaining  boundaries.  For  electron  transport,  the  electric  potential 
is  set  to  zero  at  the  land,  while  the  no-flux  condition  is  assigned  to 
all  the  remaining  boundaries.  For  ionic  charge  conservation,  the 
ionic  potential  is  set  to  0m> o  at  the  membrane/CL  interface  and  all 
the  other  boundaries  are  defined  as  no-flux.  The  parameters  for  the 
baseline  case  are  listed  in  Table  2. 

3.  Results  and  discussion 

In  Ref.  1]  we  studied  the  effect  of  CL  composition  on  perfor¬ 
mance  using  our  new  sphere-packing  approach  to  determine  the 
effective  agglomerate  surface  area  which  is  required  to  accurately 
calculate  the  diffusion  losses.  The  model  was  validated  in  Ref.  [1  by 
comparing  the  predicted  polarization  curve  against  experimental 
data.  It  was  also  shown  that  the  predicted  current  distribution 
matched  the  results  from  previously  published  numerical  studies. 
In  addition,  the  model's  predictions  for  the  effect  of  NFP  on  current 
density  at  a  constant  operating  voltage  was  shown  to  agree  with 
experimental  data.  All  of  the  results  in  Ref.  [1]  pertained  to 
spatially-uniform  CL  compositions.  Here,  we  employ  the  model  of 
[  1  ]  to  investigate  the  effect  of  graded  CL  composition  on 
performance. 

It  is  evident  that  the  concentration  and  transport  of  reacting 
species  (ions,  electrons,  O2)  in  the  CL  is  not  spatially  uniform  either 
in  the  through-thickness  or  in-plane  directions.  For  instance,  the 
proton  flux  is  the  highest  at  the  membrane/CL  interface  and  drops 
to  zero  at  the  GDL/CL  interface.  Similarly,  the  reactant  flux  is  highest 
at  the  GDL/CL  interface  and  drops  to  zero  at  the  membrane/CL 
interface.  Consequently,  the  reaction  rate  is  also  not  spatially  uni¬ 
form  within  the  CL.  These  observations  suggest  that  more  efficient 
reactions  and  higher  catalyst  utilization  could  be  achieved  by 
functionally  grading  the  CL  composition  to  fit  the  requirements 
imposed  by  these  non-uniform  reaction  rates.  Increasing  the 
effectiveness  of  the  CL  by  grading  its  composition  requires  a  good 
understanding  of  the  effects  of  the  charge  and  species  transport 
mechanisms  on  the  reaction  rate.  In  addition,  the  reaction  rate 
distribution  and  the  overpotential  in  the  CL  are  highly  dependent 
on  the  operating  regime. 

O2  mass  fraction  distributions  for  activation  (nominal  cathode 
overpotential  (NCO)  =  0.2  V),  ohmic  (NCO  =  0.5  V),  and  mass 


Table  2 

Parameter  values  used  in  the  cathode  simulation  for  the  baseline  case. 


Pressure 

P 

2.04  atm  (30  psi) 

Temperature 

T 

353  K 

Relative  humidity 

RH 

0.7 

Platinum  loading 

mPt 

0.003  kg  m  2 

Pt  to  C  weight  ratio 

Pt|C 

0.2 

Ionomer  weight  fraction 

NFP 

0.30 

Agglomerate  radius 

ragg 

3E-6  m 

Ionomer  film  thickness 

5 

145E-9  m 

Henry's  constant 

H 

31663  Pa  m3  mol-1 

Reference  exchange  current  density 

*0,ref 

0.00105  A  ur2  for  Vc  >  0.8  V 

0.06948  A  m~2  for  Vc  <  0.8  V 

Cathode  charge  transfer  coefficient 

ac 

1  for  Vc  >  0.8  V 

0.61  for  Vc  <  0.8  V 

Electron  conductivity 

0’s 

150  S  m"1 

Diffusivities 

do2,n2 

1.86E-5  m2  s”1 

dh2o,n2 

2.58E-5  m2  s'1 

Dq2,h2o 

2.47E-5  m2  s’1 

transport  (NCO  =  0.9  V)  regimes  are  shown  in  Fig.  3  for  the 
spatially-uniform  CL  composition  case.  In  the  activation  regime,  the 
O2  concentration  throughout  the  CL  is  high  and  uniform  as  the  O2 
consumption  rate  is  small  compared  to  the  rate  of  diffusion 
(Fig.  3a)).  In  contrast,  the  reaction  rate  greatly  exceeds  the  rate  of 
diffusion  in  the  mass  transport  regime  due  to  which  the  under-land 
region  experiences  greatly  reduced  O2  concentrations  (Fig.  3c)). 

The  overpotential  is  the  difference  between  electron  and  ionic 
potentials.  It  conveys  information  about  the  effects  of  both  the 
electron  and  ion  transport.  Fig.  4  illustrates  the  overpotential  dis¬ 
tribution  in  the  CL  for  activation  (NCO  =  0.2  V),  ohmic 
(NCO  =  0.5  V),  and  mass  transport  (NCO  =  0.9  V)  regimes.  The 
overall  trend  in  Fig.  4  indicates  that  the  overpotentials  are  always 
higher  under  the  land  (due  to  shorter  electron  paths)  and  close  to 
the  membrane/CL  interface  (due  to  higher  ion  concentrations).  In 
the  activation  regime,  the  spatial  differences  in  overpotential  are 
very  small;  the  overpotential  is  slightly  higher  under  the  land  than 
under  the  channel  due  to  shorter  electron  paths.  Furthermore,  the 
electrolyte  potential  only  varies  in  the  through-thickness  direction 
in  the  activation  regime  as  can  be  seen  in  Fig.  5a).  In  the  ohmic  loss 
region,  the  spatial  variations  in  overpotential  become  greater 
mainly  due  to  a  higher  ionic  potential  drop  as  illustrated  in  Fig.  5b). 
Fig.  5b)  indicates  that  the  ionic  potential  greatly  decreases  under 
the  channel  due  to  an  abundance  of  reactant  which  results  in 
higher  ion  consumption  under  the  channel.  In  the  mass  transport 
regime,  the  overpotential  distribution  in  Fig.  4c)  closely  follows  the 
ionic  potential  distribution  Fig.  5c)  which  shows  that  ions  cannot  be 
consumed  under  the  land  at  high  reaction  rates  due  to  a  lack  of 
reactant. 

The  reaction  rate  distribution  for  the  spatially-uniform  CL 
composition  is  shown  in  Fig.  6  for  the  same  three  operating  re¬ 
gimes.  The  location  of  the  highest  reaction  rate  can  be  explained  by 
combining  the  effects  of  mass  and  charge  transport  presented  in 
Figs.  3-5.  The  reaction  rate  distribution  in  Fig.  6a)  is  almost  iden¬ 
tical  to  the  overpotential  distribution  in  Fig.  4a)  which  implies  that 
the  reactions  are  dominated  by  local  overpotentials  in  the  activa¬ 
tion  region  as  the  reactant  concentration  is  almost  uniform 
throughout  the  CL.  In  the  ohmic  loss  regime,  the  highest  reaction 
rate  is  observed  at  the  junction  of  the  land  and  the  channel  at  the 
membrane/CL  interface  (Fig.  6b)).  Thus,  it  can  be  inferred  that  the 
ohmic  loss  region  is  almost  equally  dominated  by  both  reactant  and 
charge  transport  such  that  the  highest  reaction  rate  is  obtained 
where  all  species  can  be  transported  easily.  When  the  reaction  rate 
increases  further,  reactant  transport  begins  to  dominate  and  the 
highest  reaction  rate  is  obtained  under  the  channel  (Fig.  6c)).  In 
Fig.  6c),  it  is  also  observed  that  the  reaction  rate  changes  only  in  the 
in-plane  direction  under  the  land  due  to  severe  reactant  starvation 
there.  Furthermore,  close  to  the  GDL/CL  interface  in  Fig.  6,  higher 
reaction  rates  are  observed  at  the  channel/land  junction  which  can 
be  attributed  to  shorter  electron  paths. 

It  is  important  to  note  that  the  highest  reaction  rates  are  always 
obtained  at  the  membrane/CL  interface  showing  the  importance  of 
ion  concentration  on  reactions.  This  observation  suggests  that  the 
CL  composition  could  be  graded  to  improve  ion  transport. 
Accordingly,  we  first  investigate  the  case  with  graded  NFP  in  the 
through-thickness  direction  while  fixing  the  total  Nation  content. 
The  CL  is  divided  into  three  layers  of  equal  thickness,  each  with  a 
different  NFP,  but  with  a  constant  average  NFP  of  30%.  The  resulting 
performance  curves  are  presented  in  Fig.  7. 

Fig.  7  indicates  that  the  performance  is  significantly  affected  by 
the  NFP  gradient  in  through-thickness  direction.  A  higher  NFP  close 
to  the  membrane  results  in  improved  performance  as  was  experi¬ 
mentally  observed  in  Ref.  11  .  Most  of  the  improvement  is 
observed  in  the  ohmic  loss  region  because  higher  NFP  close  to  the 
membrane  improves  the  transport  of  ions  migrating  into  the  CL 
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Fig.  3.  02  mass  fraction  distributions  in  the  CL  for  different  operating  regimes  (spatially-uniform  CL  composition). 
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Fig.  4.  Absolute  values  of  overpotential  (V)  distribution  in  the  CL  for  different  operating  regimes  (spatially-uniform  CL  composition). 


Normalized  CL  width  Normalized  CL  width  Normalized  CL  width 

Fig.  5.  Ionic  potential  (V)  distribution  in  the  CL  for  different  operating  regimes  (spatially-uniform  CL  composition). 


from  the  membrane.  In  addition,  decreasing  NFP  close  to  the  GDL 
improves  reactant  transport  by  increasing  the  porosity.  Due  to 
these  effects,  the  maximum  power  density  is  improved  by  about  7% 
over  the  baseline  case.  On  the  other  hand,  if  we  increase  the  NFP 
content  close  to  the  GDL  and  decrease  the  NFP  close  to  the  mem¬ 
brane,  both  the  ion  and  reactant  transport  are  negatively  affected, 
and  hence  the  performance  decreases  as  seen  in  Fig.  7.  Further¬ 
more,  Fig.  7  shows  that  the  NFP  close  to  the  membrane  should  not 
be  increased  beyond  about  40%  due  to  diffusion  limitations.  In 
other  words,  there  is  an  optimal  NFP  gradient  in  the  through¬ 
thickness  direction.  In  Fig.  7  it  is  observed  that  diffusion  losses 
first  decrease  the  current  generation  only  in  the  mass  transport 


region  for  the  high  NFP  gradient  configuration  of  0.45_0.30_0.15, 
whereas  they  affect  the  entire  polarization  curve  for  the 
0.50_0.30_0.10  NFP  configuration. 

Next,  we  investigated  the  graded  catalyst  loading  in  through¬ 
thickness  direction.  Similar  to  the  NFP  gradient  case,  the  CL  is 
divided  into  three  layers  of  equal  thickness,  each  with  a  different 
catalyst  loading,  but  with  a  constant  average  loading  of 
0.30  mg  cm-2.  The  resulting  performance  curves  in  Fig.  8  indicate 
that  concentrating  the  catalyst  towards  the  membrane  (the  ion- 
enriched  region  in  the  CL)  increases  Pt  utilization  and  improves 
performance  in  agreement  with  the  experimental  observations  of 
[7-9  .  The  maximum  power  density  is  improved  by  about  3.2%  over 
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Fig.  6.  Reaction  rate  (mol  nr3  s)  distribution  in  the  CL  for  different  operating  regimes  (spatially-uniform  CL  composition). 


the  baseline  case.  In  contrast,  biasing  the  catalyst  concentration 
towards  the  GDL  lowers  the  performance  as  well  as  catalyst  utili¬ 
zation  as  the  reaction  rates  are  always  lower  at  the  GDL/CL  inter¬ 
face.  Attempting  to  improve  Pt  utilization  by  further  concentrating 
catalyst  at  the  membrane/CL  interface  is  not  feasible  because 
excessive  catalyst  loading  can  block  the  pores  and  reduce  reactant 
access  leading  to  diffusion  losses  as  observed  in  Fig.  8. 

We  also  studied  the  effect  of  graded  Pt|C  weight  ratio  in  the 
through-thickness  direction  while  fixing  the  Pt  loading,  NFP  and 
overall  Pt|C  ratio.  The  resulting  polarization  curves  are  illustrated  in 
Fig.  9.  It  is  important  to  note  that,  for  the  same  NFP  and  Pt  loading,  a 
higher  Pt|C  weight  ratio  implies  lower  Nation  and  carbon  loadings 
which  result  in  higher  porosity.  Hence,  higher  Pt|C  would  improve 
reactant  transport,  whereas  lower  Pt|C  would  benefit  the  ohmic 
region,  -ig.  9  shows  that  a  lower  Pt|C  close  to  the  membrane 


Fig.  7.  Effect  of  NFP  grading  in  the  through-thickness  direction  on  cathode  perfor¬ 
mance;  the  average  NFP  is  fixed  at  30%. 


Fig.  8.  Effect  of  catalyst  grading  in  the  through-thickness  direction  on  cathode  per¬ 
formance;  the  average  catalyst  loading  is  fixed  at  0.30  mg  cm-2. 


improves  performance  in  agreement  with  our  previous  findings 
because  a  lower  Pt|C  corresponds  to  a  higher  Nafion  loading  which 
improves  ion  transport.  Concurrently,  a  higher  Pt|C  close  to  the  GDL 
enhances  porosity  and  improves  diffusion.  The  maximum  power 
density  is  improved  by  about  12%  over  baseline  for  the 
0.15_0.25_0.35  Pt|C  configuration.  This  performance  improvement 
is  higher  than  that  seen  in  Fig.  7  due  to  the  combined  effects  of 
highly  improved  reactant  transport  (higher  porosity)  close  to  the 
GDL  and  higher  Nafion  loading  close  to  the  membrane.  On  the  other 
hand,  the  inverse  configuration  lowers  performance  as  shown  in 
Fig.  9. 

In  addition,  we  studied  the  case  for  varying  Pt|C  ratio  but  with 
fixed  Pt  and  Nafion  loadings.  In  this  case,  varying  Pt|C  only  changes 
the  carbon  loading  and  the  porosity.  The  performance  curves  in 
Fig.  10  indicate  that  the  diffusion  loss  region  is  improved  for  both 
configurations.  It  should  be  noted  that  the  average  porosity  de¬ 
creases  for  both  configurations  from  0.55  to  around  0.5  which  in¬ 
creases  the  effective  surface  area  and  improves  performance  in  the 
diffusion  loss  region.  The  ohmic  loss  region  is  not  noticeably 
improved  for  either  configuration  as  the  electron  potential  drop  in 
the  through-thickness  direction  is  insignificant.  Fig.  10  indicates 
that  a  greater  performance  improvement  is  observed  for  the 
configuration  with  higher  Pt|C  ratio  close  to  the  GDL;  higher  Pt|C 
close  to  the  GDL  provides  more  pathways  for  reactant  access  to  the 
active  area  and  improves  mass  transport. 

Another  important  observation  from  Figs.  3-6  is  that  the  reac¬ 
tion  rates  are  always  higher  under  the  channel  than  under  the  land 
in  the  mass  transport  regime.  This  suggests  that  the  CL  composition 
could  be  functionally  graded  in  the  in-plane  direction  to  improve 
reactant  transport. 


Fig.  9.  Effect  of  Pt|C  grading  in  the  through-thickness  direction  on  cathode  perfor¬ 
mance  for  constant  catalyst  loading  and  NFP.  The  average  Pt|C  is  fixed  at  0.25. 
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Fig.  10.  Effect  of  Pt|C  grading  in  the  through-thickness  direction  on  cathode  perfor¬ 
mance  for  constant  catalyst  and  Nafion  loadings  (mPt  =  0.3  mg  cur2, 
mN  =  0.514  mg  cm-2).  The  average  Pt|C  is  fixed  at  0.25. 
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Fig.  12.  Effect  of  catalyst  grading  in  the  in-plane  direction  on  cathode  performance; 
the  average  catalyst  loading  is  fixed  at  0.30  mg  cm  2. 


We  investigate  this  idea  by  first  grading  the  NFP  in  the  in-plane 
direction  while  fixing  the  average  NFP  in  the  CL  at  30%.  We  consider 
the  case  where  the  CL  is  divided  into  two  parts  (under-channel  and 
under-land)  each  with  different  NFP  values.  The  resulting  polari¬ 
zation  curves  in  Fig.  11  show  that  higher  NFP  under  the  channel 
improves  performance  in  the  diffusion-loss  region.  Increasing  the 
under-channel  NFP  improves  ion  conductivity  in  the  through¬ 
thickness  direction  leading  to  higher  ion  transport  towards  the 
GDL  interface  under  the  channel  where  the  reactant  concentration 
is  high  enough  to  consume  these  ions.  Increasing  the  under¬ 
channel  NFP  improves  the  power  density  by  about  6.5%  over 
baseline.  In  contrast,  increasing  the  NFP  under  the  land  results  in 
poor  performance  especially  in  the  diffusion  loss  region  because 
the  under-land  region  suffers  from  insufficient  O2  concentration  to 
consume  the  ions.  Fig.  11  also  indicates  that  increasing  the  under¬ 
channel  NFP  beyond  about  40%  exacerbates  diffusion  losses.  The 
effect  of  decreased  porosity  is  first  felt  only  in  the  diffusion  loss 
region  for  NFP_0.45_0.15  configuration;  further  NFP  grading  causes 
diffusion  losses  throughout  the  polarization  curve  as  see  for  the 
NFP_0.50_0.10  configuration. 

Next,  we  examine  the  effect  of  catalyst  grading  in  the  in-plane 
direction  to  improve  reactant  transport  in  the  mass  transport 
regime.  Similar  to  the  laterally-graded  NFP  case,  we  divided  the  CL 
into  under-channel  and  under-land  regions,  each  with  a  different 
catalyst  loading.  The  effect  of  in-plane  catalyst  gradient  on  perfor¬ 
mance  in  Fig.  12  shows  that  concentrating  the  catalyst  under  the 
channel  improves  performance  in  the  diffusion-loss  region.  The 
improvement  can  be  attributed  to  increased  catalyst  utilization 
because  the  under-channel  region  is  enriched  with  reactant 
compared  to  the  under-land  region  especially  in  the  mass  transport 
regime.  The  maximum  power  density  is  increased  about  2.5%  over 
baseline  by  concentrating  the  catalyst  under  the  channel.  Similar  to 


previous  results,  there  is  an  optimum  gradient  for  the  Pt  loading 
due  to  diffusion  limitations. 

Next,  we  study  the  effect  of  Pt|C  grading  in  the  in-plane  direc¬ 
tion  while  keeping  the  Pt  loading  and  NFP  content  constant  in  the 
CL.  Again,  the  under-land  and  under-channel  regions  of  the  CL  are 
loaded  with  different  Pt|C  ratios.  The  resulting  polarization  curves 
presented  in  -ig.  13  indicate  that  a  lower  Pt|C  ratio  under  the 
channel  improves  performance  especially  in  the  diffusion-loss  re¬ 
gion.  The  performance  improvement  in  Fig.  13  can  be  attributed  to 
higher  Nafion  and  carbon  loadings  under  the  channel  corre¬ 
sponding  to  a  lower  Pt|C  ratio;  then,  improved  ion  conduction 
under  the  channel  increases  performance.  The  maximum  power 
density  is  improved  about  9.5%  over  baseline  as  shown  in  Fig.  13. 
Although  Fig.  13  reveals  the  same  trends  as  in  Fig.  11,  the  former 
yields  higher  performance  improvement  due  to  better  reactant 
transport  under  higher  porosity.  As  expected,  the  inverse  configu¬ 
ration  results  in  lower  current  densities  in  the  mass  transport 
regime. 

In  ig.  13  the  effect  of  graded  Pt|C  is  dominated  by  the  effect  of 
Nafion  loading  (similar  to  the  case  in  Fig.  11 )  which  varies  due  to  the 
constant  NFP  assumption.  Therefore,  we  next  isolate  the  effect  of 
laterally  grading  the  carbon  loading  by  varying  the  Pt|C  ratio  at  a 
fixed  Nafion  loading.  The  resulting  performance  curves  presented 
in  Fig.  14  show  that  the  configuration  with  higher  Pt|C  under  the 
land  does  not  experience  a  noticeable  mass  transport  loss  due  to  a 
locally  lower  carbon  loading  and  higher  porosity.  Hence,  such  a  Pt|C 
grading  improves  reactant  transport  under  the  land  which  typically 
suffers  from  diffusion  losses  at  high  reaction  rates.  In  contrast,  the 
inverse  case  of  lower  Pt|C  under  the  land  reveals  a  distinct  mass 
transport  loss  region.  However,  the  case  of  lower  Pt|C  under  the 
land  performs  slightly  better  in  the  ohmic  loss  region  as  the 
effective  under-land  surface  area  is  increased  over  the  baseline 
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Fig.  11.  Effect  of  NFP  grading  in  the  in-plane  direction  on  cathode  performance;  the 
average  NFP  is  fixed  at  30%. 


Fig.  13.  Effect  of  Pt|C  grading  in  the  in-plane  direction  on  cathode  performance  for 
constant  catalyst  loading  and  NFP.  The  average  Pt|C  is  fixed  at  0.25. 


F.C.  Cetinbas  et  al.  /  Journal  of  Power  Sources  270  (2014)  594-602 


601 


Fig.  14.  Effect  of  Pt|C  grading  in  the  in-plane  direction  on  cathode  performance  for 
constant  catalyst  and  Nafion  loadings  (mPt  =  0.3  mg  crrr2,  mN  =  0.514  mg  cm  2).  The 
average  Pt|C  is  fixed  at  0.25. 

case.  Higher  effective  surface  area  under  the  land  increases  current 
generation  in  the  ohmic  loss  regime  leading  to  higher  reaction  rates 
under  the  land  as  illustrated  in  Fig.  6b).  Due  to  this  improvement  in 
the  ohmic  loss  region,  the  maximum  power  density  is  increased  by 
about  4%  over  baseline. 

So  far,  we  studied  the  effects  of  graded  CL  composition  in  the 
through-thickness  and  in-plane  directions  separately.  Graded 
composition  in  the  through-thickness  direction  mainly  improves 
performance  in  the  ohmic  loss  regime  by  improving  ion  transport. 
Graded  CL  composition  in  the  in-plane  direction  mainly  improves 
performance  in  the  mass  transport  regime  by  improving  reactant 
transport.  It  is  worthwhile  to  investigate  if  grading  the  CL  compo¬ 
sition  in  both  directions  simultaneously  can  provide  a  performance 
that  exceeds  the  benefits  of  CL  grading  solely  in  either  individual 
direction. 

Accordingly,  we  first  examined  the  benefit  of  grading  NFP  in  the 
through-thickness  direction,  the  in-plane  direction,  and  bidirec¬ 
tionally  as  illustrated  in  Fig.  15  while  keeping  the  average  NFP 
constant  at  30%.  Fig.  15  indicates  that  grading  the  NFP  bidirec¬ 
tionally  improves  performance  both  in  the  ohmic  and  mass  trans¬ 
port  loss  regimes.  Bidirectional  grading  improves  ion  transport  in 
the  through-thickness  direction  concentrating  Nafion  close  to  the 
membrane.  At  the  same  time,  concentrating  Nafion  under  the 
channel  improves  ion  transport  and  reactant  consumption  in  the 
mass  transport  regime  for  which  the  reactant  is  abundant  only 
under  the  channel. 

Finally,  we  compare  the  performance  obtained  by  grading  the 
catalyst  in  the  through-thickness  direction,  the  in-plane  direction, 
and  bidirectionally  as  illustrated  in  Fig.  16  while  keeping  the 
average  catalyst  loading  constant  at  0.30  mg  cm-2.  As  expected,  the 
bidirectional  grading  superposes  the  positive  effects  of  the 


Fig.  15.  Effect  of  bidirectional  NFP  grading  on  performance;  the  average  NFP  is  fixed  at 
30%. 
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Fig.  16.  Effect  of  bidirectional  catalyst  grading  on  performance;  the  average  catalyst 
loading  is  fixed  at  0.30  mg  cm  2. 


through-thickness  and  in-plane  cases.  Therefore,  bidirectionally- 
graded  catalyst  loading  provides  improved  performance  in  both 
ohmic  and  diffusion-loss  regions. 

Figs.  15  and  16  indicate  that  bidirectional  grading  would  be  a 
desirable  option  if  such  a  CL  could  be  manufactured.  In  this  study, 
we  have  only  considered  one  out  of  several  possible  multi¬ 
directional  grading  configurations.  The  results  are  promising,  and 
suggest  that  performance  could  be  improved  further  by  simulta¬ 
neously  optimizing  the  CL  composition  gradients  along  multiple 
directions.  Further  study  is  required  to  optimize  the  CL  topology  for 
superior  performance  over  a  range  of  operating  conditions. 

4.  Summary  and  conclusions 

In  this  study,  we  examined  functionally-graded  catalyst  layer 
compositions  by  employing  our  two-dimensional  cathode  model 
which  incorporates  an  improved  agglomerate  approach  [1]  to 
accurately  model  diffusion  limitations.  First,  the  effects  of  charge 
and  mass  transport  mechanisms  on  reaction  rate  distribution  were 
investigated  for  different  operating  regimes.  It  was  observed  that 
the  highest  reaction  rates  are  always  obtained  at  the  membrane/CL 
interface  for  all  operating  regimes  due  to  higher  local  ion  concen¬ 
trations.  In  addition,  higher  reaction  rates  were  observed  under  the 
channel  in  the  mass  transport  regime  due  to  the  greater  abundance 
of  reactant. 

Based  on  these  observations,  we  studied  graded  CL  composi¬ 
tions  to  improve  ion  and  reactant  transport  in  the  CL.  A  higher  NFP 
close  to  the  membrane  improves  ion  transport  and  performance  in 
the  ohmic  loss  region.  Higher  NFP  under  the  channel  improves 
performance  in  the  mass  transport  region  due  to  both  improved  ion 
and  reactant  transport.  Concentrating  Pt  close  to  the  membrane 
improves  performance  by  increasing  catalyst  utilization  in  this  ion- 
rich  region.  In  addition,  higher  catalyst  loading  under  the  channel 
improves  performance  in  the  mass  transport  regime  in  this 
reactant-rich  region.  Optimum  gradients  were  found  for  both  NFP 
and  catalyst  loading.  Increasing  the  gradient  beyond  a  critical  value 
results  in  decreased  performance  due  to  diffusion  limitations. 

We  also  studied  bidirectionally-graded  CL  compositions  for  the 
first  time.  Bidirectional  grading  yields  superior  performance  over 
unidirectional  grading  due  to  simultaneous  improvements  in  both 
the  ohmic  loss  and  mass  transport  regimes. 
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Nomenclature 

aagg  effective  surface  area  in  CL  m2  m  3 

Co  ionomer  film  surface  oxygen  concentration  mol  nrT3 

Cs  agglomerate  surface  oxygen  concentration  mol  m”3 

Dij  binary  diffusivities  of  i  and  j  m2  s”1 

D?Jf  effective  binary  diffusivities  of  i  and  j  m2  s”1 

D  diffusivity  of  oxygen  in  ionomer  m2  s”1 

Deff  effective  diffusivity  of  oxygen  in  ionomer  m2  s”1 

Dp  particle  diameter 

F  Faraday's  constant  C  mol”1 

/  overlapping  ratio 

H  Henry's  constant  atm  m3  mol-1 

I„L  current  generation  per  unit  volume  of  CL  A  m”3 

i^ef  reference  exchange  current  density  A  cm-2 

Mn  molecular  weight  of  gas  mixture  kg  mol”1 

Mj  molecular  weight  of  species  j  kg  mol”1 

m  number  of  contact  points 

mpt  Pt  loading  mg  cm”2 

NFP  Nation  weight  ratio 

n  number  of  agglomerate  per  unit  volume  of  CL 

P  pressure  atm 

P0l  oxygen  partial  pressure  atm 

Ph2  hydrogen  partial  pressure  atm 

Pt|C  platinum  carbon  weight  ratio 

R  gas  constant  mol  cm”2  s”1 

Rq2  volumetric  oxygen  reaction  rate  mol  m”3  s”1 

PHo2  volumetric  reaction  rate  for  water  vapor  mol  m”3  s”1 

ragg  agglomerate  radius  m 

rPt  Pt  particle  radius  m 

T  temperature  I< 

Th  Thiele  modulus 

tci  catalyst  layer  thickness  m 

ac  charge  transfer  coefficient 

am  packing  factor 


aoD  electroosmotic  drag  coefficient 
5  ionomer  film  thickness  m 

cgdl  porosity  of  GDL 

ev,cr  critical  porosity 

ev,o  limiting  porosity 

£v  porosity  of  CL 

6-agg  ionomer  volume  fraction  in  agglomerate 

p  density  kg  m”3 

a s  electron  conductivity  S  m”1 

effective  electron  conductivity  S  m”1 
om  ionic  conductivity  S  m”1 

effective  ionic  conductivity  S  m”1 
{ f>s  electron  potential  V 

(j)m  ionic  potential  V 

cj)RU  local  relative  humidity 

oji  mass  fraction  of  species  i 
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